High Energy Dynamics of the Single Impurity Anderson Model 
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The quantitative control of the dynamic correlations of single impurity Anderson models is es- 
sential in several very active fields. We analyze the one-particle Green function with a constant 
energy resolution by dynamic density-matrix renormalization. In contrast to other approaches, 
sharp dominant resonances at high energies are found. Their origin and importance are discussed. 
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Single impurity models are at the very basis of the de- 
scription of strong correlation phenomena. Landmarks 
are the Kondo problem! and the single impurity Ander- 
son model (SIAM)£, for a review see Ref. 0. 

The interest in the quantitative analysis of SIAMs has 
been intensified considerably by the advent of a system- 
atic mapping of strongly correlated lattice models to ef- 
fective SIAMs supplemented by a self-consistency condi- 
tion. This is the key point of dynamic mean-field theory 
(DMFT)2£ which is based on an appropriate scaling of 
the non-local parts of the Hamiltonian&i, for reviews see 
Ref. |33- In recent years, the DMFT is applied very suc- 
cessfully in combination with ab-initio density functional 
calculationsiSiii. In this way, the unbiased knowledge 
about the bands could be enhanced by the inclusion of 
interaction effects between the excited quasi-particles. It 
turned out that the combination of density functional re- 
sults and DMFT makes the quantitative understanding 
of spectroscopic data possible^. 

So far, the methods applied to the SIAM were de- 
signed to capture the low-energy physics, in particular 
the fixed points and the thermodynamicsi&il. The nu- 
merical renormalization group (NRG) was later extended 
to calculate also dynamic, i.e., spectral information. It 
provides reliable data on the scale of the Kondo tempera- 
tures Tk , see Ref. I3lll5l and references therein. On larger 
scales, the energy resolution is less well-controlled. 

But in various applications the behavior at higher en- 
ergies is important to achieve quantitative accuracy. For 
instance, the self-consistency cycle of the DMFT mixes 
modes at all energies. Hence, excellent quantitative con- 
trol over the dynamics at high energies is indispensable, 
even if finally only the behavior at low energies matters. 

Another application is the optical control of isolated 
5=1/2 impurities or quantum dots coupled to narrow 
bands. If the impurities differ so that the energy between 
the singly occupied ground state and the excited double 
occupancy differs, they can be switched selectively from 
the ground state to the double occupancy (and back) by 
shining light at the resonant frequency onto the sample. 
The life-time of the double occupancy, i.e., the inverse 
line width of the resonance, determines how well the res- 
onance condition has to be met, how selective the switch- 
ing can be, and how stable the excited state is. 

In view of the above, we perform a numerical investi- 



gation which aims to describe both the low-energy dy- 
namics and the high-energy dynamics quantitatively. To 
this end, we use an energy resolution which is constant 
for all energies. Features at low energies are not as deli- 
cately resolved as by NRG, but in return features at high 
energies are much better under control. We apply the dy- 
namic density-matrix renormalization (D-DMRG)i£*iLi£ 
to compute the one-particle propagator. The DMRG is 
a real-space approacbi^iSSiSi which works best for open 
boundary conditions so that it is particularly well-suited 
to treat impurity problems. 

The model investigated at zero temperature is the sym- 
metric Anderson model 

H = ^ f-dn d ,a + U n d ^n d ^ + V ^(<4c 0(T + n - c -) 

<7 <7 
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with arbitrary density of states (DOS) po(u>) of the U = 
one-particle Green function Gq(u>) of the d-electron. The 
parametrization in Q is chosen such that the coefficients 
(£nj7n) are the continued fraction coefficients of the hy- 
bridization function 
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This model has particle- hole symmetry iff = — U/2 and 
e n = for all N > 0. The representation as continued 
fraction 2 ^ (see Fig. QJi) is optimum for the DMRG which 
is designed for chains. We look at a generic situation with 
finite band width W = 2D. For simplicity we choose a 
T(w) with semi-elliptic DOS, i.e., 7„ = D/2. For V = 
D/2, the free DOS p Q (uj) = 2^D 2 - uj 2 /(itD 2 ) is also 
semi-elliptic. 

The problem illustrated in Fig. is mapped by two 
standard Jordan- Wigner transformations to two XY spin 
1/2 chains, the S'-chain and the T-chain. The S-chain 
results from the f fermions, the T-chain from the J, 
fcrmions. They are coupled at site zero where the 
density-density coupling is mapped to the product of 
z-components. The resulting chain is depicted for the 
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FIG. 1: (a) Single impurity model with the bath as half- 
infinite chain, (b) equivalent spin model: T-spins come from 
f -fermions, S'-spins from j-fermions. 



symmetric SIAM in Fig. \T]p. The couplings are given by 
^0 = to' = V and U = = 7; for i > 1. The map- 
ping from fermions to spins avoids the fcrmionic Fock 
space which would imply numerically difficult long-range 
effects. The mapping makes the Hilbert space the direct 
product of the local Hilbert spaces at each site. 

The DMRG can easily determine the ground state |0) 
and its energy Eq for a finite chain. So the chain in 
Fig. \T]p is truncated such that there are L spins in the 
upper and in the lower part of the chain corresponding 
originally to a truncated bath of L — 1 fermions plus the 
impurity. The dynamic quantity we are interested in is 
the retarded Green function at zero temperature 

G>(uj + irj) = (Q|So (w + ir)-(H- i^r'S^O) (3) 

where the superscript > implies that represents only 
the part of the usual Green function at non-negative fre- 
quencies. In the symmetric case, the complete function is 
recovered by G(u> + irj) = G > (u> + irj) — G > (— uj — irj). In 
the asymmetric case, G < (uj + irj) = (0|S + (w + irj— {H — 
Eq))~ 1 Sq |0) must be determined separately whereby 
G(uj + irj) = G > (uj + irj) — G < (—uj — irj) is obtained. We 
stress that G(uj + irj) is the fermionic propagator even 
though it is computed in terms of spins after the Jordan- 
Wigner mapping. 

The key idea of the dynamic DMRG is to include the 
real and the imaginary part of the correction vector |c) 
in the target states of a standard DMRG algorith: 
The natural choice is |c) = (uj + if] — (H — E )) 
The computation of |c) is numerically the most demand 
ing step due to the inversion of an almost singular non- 
hcrmitean matrix. We prefer to stabilize this inversion by 
optimized algorithms^ instead of using the variational 
approach proposed by Jeckelmann 2 4 which requires a 
minimization in a high-dimensional Hilbert space. 

The numerical calculations cannot be performed for 
r\ = 0. Even small values of 77 are very time consuming. 
So we compute first G(uj + irj) at finite rj. The spectral 
density p^{uj) = —^QG(uj + irj) can be seen as the ac- 
tual spectral density p(oj) convoluted by the Lorentzian 
p L (w) = (rj / tt) / (uj 2 +?7 2 ) of width ?/. Hence it is possible 
to retrieve p(w) by deconvolution. A standard technique 
for deconvolution is Fourier transformation, realized best 
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by fast Fourier transforms, division by exp(— r/r) plus 
low-pass filtering, and the inverse transform. A flexible 
alternative with similar properties is the explicit matrix 
inversion of the convolution procedure 2 ^. 




FIG. 2: (a) Spectral densities for V = D/2 broadened by 
T] = 0.1D. Chain length L = 80 fermionic sites; 128 DMRG 
states kept, (b) Spectral densities from (a) deconvoluted in 
the time domain. NRG data (A — 1.8, 1500 states) depicted 
by thin dashed lines. The exact curves represent the analytic 
L — <x results. 



In Fig. generic broadened spectral densities are 
plotted as they are computed by D-DMRG. Obviously, 
the value p^(0) is not independent of U. Increasing the 
chain length L does not lead to any significant change in 
the data (not shown) . Fig. [2b displays the deconvoluted 
data. The deconvolution works very well except for some 
slight overshooting in regions where the spectral density 
changes rather abruptly. In particular, the value p(0) 
is pinned to D/(2ttV 2 ) independent of U as required by 
Friedel's sum rule and the density of states rule^SSiSL 2 ^. 
We take this fact as convincing evidence for the reliability 
of the numerical algorithm. 

The central peak at u = is the Abrikosov-Suhl 
resonance (ASR). For larger U (smaller V) its width 
decreases rapidly so that the ASR is very difficult to 
resolve^ unless more elaborate deconvolution schemes 
are devised^. So a quantitative analysis of the ASR is 
postponed to future work. 

For comparison, the thin dashed lines in Fig.^ depict 
standard NRG datai^i. For small frequencies NRG is 
well-controlled. Indeed, for \oj\ < D/3, NRG and D- 
DMRG data agree excellently lending further support to 
the D-DMRG approach. Outside the band, the NRG 
spectra appear to be too wide due to the chosen constant 
broadening on a logarithmic mesh. This broadening does 
not account for the absence of states outside the bare 
band. The NRG does not possess intrinsic information 
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about the peak widths. The position of the high energy 
peak in the raw NRG data, however, coincides with the 
D-DMRG result. 

An increase in U leads to the formation of Hubbard 
satellites below and above the free band (Fig. 0). They 
are situated at energies w up /iow = ±(£//2 + <5 s hift), £>shift > 
and become more pronounced on increasing U in two 
ways. They capture more weight and they become 
sharper. For Tk — > the weight reaches 1/2, see Ref.1281 
The sharpening has not been discussed quantitatively be- 
fore although the extended non-crossing approximation 32 
provides sharp satellites if they lie outside the bare bands, 
see e.g. Fig. 1 in Ref. |2^. Recently, indications have 
occured 2 ^ that other standard algorithms overestimate 
the width of the Hubbard satellites. The exaggerated 
width of the NRG data at high energies results from the 
Gaussian broadening of the order of the energy ranged. 

To investigate the line shapes of the satellites we plot 
them for various values of U and V in Fig. [21 The ASR 
at lo = is not displayed since it is too much smeared 
out at r\ = 0.1D for larger values of the interaction. The 
shifts Jshift increase on increasing V; they decrease on 
growing interaction U. The widths behave qualitatively 
similar. A complete deconvolution suffers unfortunately 
from severe overshooting due to the sharpness of the res- 
onance. To make the analysis nonetheless quantitative 
we fit the broadened data by Lorentzians plus an offset 
(Fig. 0) • These fits work very well for large values of U 
and not too large values of V. 




FIG. 3: D-DMRG data of the upper Hubbard satellite at 
77 = O.ID. Thin lines are fits with Lorentzians and an offset 
oc ?7cff/(^cff + ( w — w up ) 2 ) + C. The fits were done for the 
intervals shown, (a) dependence on U at constant V (L = 80); 
(b) dependence on V at constant U (L = 40). 
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FIG. 4: Widths (symbols) of the Hubbard satellites as found 
from the fits in Fig. [3] The fit intervals are given in the legend. 
The lines show various power law fits, (a) dependence on U 
at V = D/2; (b) dependence on V. 



width rjes of the convolution of two Lorentzians of widths 
771 and 772 is r] e g = r\\ + r\%. From the effective widths 
?7cff we deduce the true half-widths at half-maximum 
(HWHM) of the Hubbard satellite by subtracting the ar- 
tificial broadening 77, i.e., HWHM = rj e gf — rj. In Fig. 21 
the widths are depicted as function of U and of V. The 
results show that the HWHM are proportional to V 4 /U 2 . 
The deviations for smaller widths must be attributed to 
the numerical constraints, e.g., finite 77 and finite chain 
length L. The deviations for larger widths, mainly for 
larger values of V and smaller values of U result from 
the vicinity of the bare bands. Fig. displays the analo- 
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To deduce the true width of the Hubbard satellite we 
assume that it is well described by a Lorentzian. The 



FIG. 5: Position shifts of the Hubbard satellites analyzed 
like the widths in Fig.0 (a) V = D/2; (b) see legend. 
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gous analysis for the shifts S s hiit of the satellite positions. 
Again, strong evidence for power law behavior is found, 
namely S shiit (xV 2 /U. 

How can the above findings be understood? Let us 
start by the positions. The energy levels of isolated im- 
purities, i.e., V — are at ±E//22A. Switching on V mixes 
the impurity levels with the bath states which lie in the 
interval [—D,D]. If U is large compared to D second 
order perturbation in V implies that the impurity levels 
are repelled from the bath states. The shift should be of 
the order of J = 8V 2 /U, sec Eq. (11) in Ref. 34, which 
agrees nicely with the power laws in Fig. [3J 

The widths of the satellites have been discussed quan- 
titatively when they lie within the bare band^. If the 
satellites lie outside, perturbation theory in U, to second 
order or the random phase approximation, implies that 
a finite width is to be expected at least for U < 6D. But 
the reasoning in powers of U is not helpful for U > 6D 
and it does not explain the power laws found. So we 
return to perturbing in powers of V. The impurity lev- 
els mix with particle-hole excitations in the bands, see 
Eq. (10a) in Ref. |3 In the symmetric case the dou- 
bly occupied electron and hole state are degenerate so 
that mixing with particle-particle (or hole-hole) states 
matters also, see Eq. (12) in Ref. yj- The mixing is of 
order J = 8V 2 /U. So Fermi's golden rule implies a life 
time of J 2 No where No measures the density of states 
with which the impurity level mixes; -/Vo is of the order 
of D -1 . Indeed, HWHM cx J 2 explains conclusively the 



data of Fig. H 

So far, the width of the Hubbard satellites for U > 2D 
was extracted under the assumption that the satellites 
are Lorentzians. Further investigations of the line shape 
are urgently called for. Numerically, improvements of 
the resolution are necessary to determine the line shape 
of the satellites explicitly. Analytically, the quantitative 
argument for the widths must be supplemented by an 
explicit calculation of the line shapes for V 2 /U — > 0. 

In summary, we have investigated the dynamic prop- 
agator of the SIAM by D-DMRG. This powerful large- 
scale algorithm provides information with a constant en- 
ergy resolution. Up to moderate interactions U s=s 2D, 
deconvolution yields the explicit spectral densities. For 
larger interactions, the width of sharp resonances can be 
extracted by fitting Lorentzians. In particular, we ana- 
lyzed the positions and widths of the Hubbard satellites. 
The shifts are of order V 2 /U due to level repulsion; the 
line widths are of order V 4 /U 2 . 

Especially the sharpness of Hubbard peaks is missed by 
other zero temperature algorithms for the SIAM. Hence 
the D-DMRG is a very valuable complementary tool. Po- 
sition and width of the Hubbard satellites are impor- 
tant for several widely used applications, e.g., in the self- 
consistency cycle of the DMFT. 
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